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1 Introduction 

In Regge theory the high-energy asymptotics of the hadron-hadron scattering aniphtudes is 
determined by singularities of partial waves in the complex angular momentum plane. It was 
observed ffl that the regularity of quantum mechanics, which relates high-energy scattering 
amplitudes to the singularities of the partial waves in the complex angular momentum plane is 
valid for quantum field theory as well. Namely, in Regge kinematics 

sy^-tr^M^ (1) 

where M is a characteristic hadronic mass scale, the hadron-hadronic scattering amplitude 

A{s,t) 

e+ioo , 

e—ioo 

is governed by singularities of the partial waves /(J, t), by Regge poles and Regge cuts. Among 
the Regge singularities there is one with vacuum quantum numbers, the so-called Pomeron, 
which provides the dominant contribution to the scattering amplitude. The bootstrap conjecture 
has been proposed, according to which all particle-like excitations correspond to some Regge 
singularity and are related to each other via unitarity of the S-matrix and sum rules. 

However the programme of building the axiomatic quantum field theory from assumptions of 
only unitarity and analyticity of the S-matrix failed, because Regge theory itself does not allow 
to calculate the positions of these singularities. Now, QCD as a theory of strong interaction is 
called to describe the Regge behavior of the scattering amplitudes p. V.N. Gribov proposed 
the conjecture, that reggeons form new collective excitations and QCD in the high-energy limit 
can be replaced by an effective reggeon theory y]. This has been confirmed in series of works, 
initiated by L.N.Lipatov ^. It was shown that in leading logarithmic approximation (LLA), 
which is the natural approximation in the Regge limit of QCD, A{s,t) can be expressed as a 
sum of Feynman diagrams describing the multiple exchange of reggeized gluons in the f-channel. 
The perturbative expressions for corresponding Feynman diagrams including large logarithmic 
factors ol"'log^s (m=n,n-l,...) have to be resummed to all orders in a^, because bare gluons 
and quarks are not a good approximation in the Regge limit. The leading contribution (m = n) 
comes from ladder diagrams, corresponding to exchange of n-reggeons in t-channel. Being built 
from an infinite number of perturbative gluons, the reggeons carry the quantum numbers of the 
gluon and become a new collective excitation in the Regge limit (1.1). It is well known that 
the leading logarithmic approximation results in an asymptotics, which violates the Froissart 
bound. Unitarity is restored by taking into account sub-leading contributions as well. In the 
generalised leading logarithmic approximation (GLLA) some minimal set of non-leading terms 
is included to restore unitarity [0]. The interaction of the reggeons is determined by LLA. 

In LLA the dominant contribution to the partonic scattering amplitude comes from the 
soft gluons and this leads to the gluon reggezation property. It can be shown that infrared 
divergences are cancelled for colourless external states due gauge invariance. 

Of course the distribution of partons inside the hadron are described by the nonperturbative 
wave function of the hadron. The non-perturbative effects can be taken into account in the 
approach of constructing the high-energy effective action Q. However, the perturbative inves- 
tigation of high-energy QCD makes sense as a first approximation to begin with. In particular 
results obtained in the BFKL pomeron approximation are in a good agreement with experimen- 
tal data of semi-hard processes and especially of deep-inelastic scattering at small x. Then, some 
information about non-perturbative corrections can be extracted from analyzing the behaviour 
of the perturbative series in the infrared region W\ . 



2 High-energy QCD as an integrable model 

As mentioned above, the high-energy asyniptotics of the scattering amphtude in leading loga- 
rithmic approximation is determined by the contribution of diagrams, describing two reggeon 
exchange in the t-channel. It is a result of summing up an infinite number of the famous ladder 
diagrams Q and corresponds to the exchange of the Pomeron. Contributions of diagrams with 
three (Odderon) and more reggeized gluons can be considered as higher corrections. Pomeron 
contribution corresponds to elastical scattering. Diagrams with Odderon exchange describe 
processes with the exchange of negative charge parity. 

The dominating contribution in LLA comes from the multi-Regge kinematics: 

s = {pA+ Psf ~ 2pA • PB (3) 

i = l,..n + l, ko=pA, ki = qi+i-qi, K+i =Pb 

s » Sj >[ qi [^ 

n 
S1S2 ■ ■ ■ Sn+1 = sW{-k\^) 
i=l 

fel • PA < A:2 • PA < ••• <^kn-pA 
kl ■ PB ^ k2 ■ PB ^ ■■■ ^ K-PB 

where k± is defined by the Sudakov decomposition 

, u k ■ PA a , k-B It , J a 

k^ = pT, H p^ + kK 

PA -PB PA- PB 

Owing to this decomposition the scattering amplitudes in the Regge limit exibit the remarkable 
separation of the longitudinal and transverse directions with respect to the plane, spanned by 
the momenta of the initial particles. 

In the generalized leading logarithmic approximation the interaction between reggeons is 
elastic and pairwise. We shall restrict ourselves to the case where the number of reggeons in 
the t-channel A^ is conserved. For a given A^ the reggeon Green function fuf.} satisfies the 
Bethe-Salpeter-like equation 

-/{u = /{2} + E41o.}A..} (4) 

Kj 

The set {z^} = {ii,..ir) labels the reggeons, i = G stands for gluon and i = F 01 i = F 
stand for reggeized quarks of corresponding helicity. The partial wave / carries also an index 
ai, labelling the gauge group representation of the corresponding reggeons and depends in their 
transverse momentum k±i or their impact parameter Xj. 

The r-reggeon contribution to the partial wave is obtained by contracting the r-reggeon 
Green function with the parton distribution functions of the scattered particles. The angular 
momentum is 

j = l + ^_!I (5) 

where rj is the number of exchanged fermions. The pairwise interaction of the reggeons is 
described by the Hamiltonian 

W^''^Vn,..v.v = j^ J dk[dk'j5{h + kj - k'j - k'j)[{T, ® r,)7^7ii,,/,,,..,,..j- ..^ + (6) 



The first term in the square bracket corresponds to the interaction via an s-channel gluon: 

m ® r.)H,„,„.,,j = n <,»;,iTna^an{Tna,a,, (7) 

and the second one corresponds to the interaction via an s-channel fermion 

m ^ T,)g^^^^^,^^ = n '^a„a,.(T-'')«,:«(r°0«a/, (8) 

with the sum over the fermion colour states a. The overall group state in the t-channel has to 
be the gauge singlet. 

While the longitudinal part of the scattering amplitude is extracted as a kinematical factor, 
the transverse dynamical part can be described by simple Feynmann rules, corresponding to 
the multi-regge effective action [|5|. These graphical rules allow the simple derivation of the 
interaction kernels, entering the equation (pi). As operators in impact space these kernels take 
the following form 

'Hgg = Hg + Hq (9) 

Hff = Hg + Hp 

Hfg = Hg + Vi2HpVi2 

Gfg = {xl2dl)-^ 

where 

Hg = -2^(1) + d^^ logxia^i + d:^^ logxia^a + logai^a (10) 

H^F = -2v(i) + ai"'+"/' iogxi2ai""/' + d^"^'^ iogxi2a^/' + log^i^s 

and "0(1) = ~1 is the Euler number. 7^i2 represents the operator permuting the reggeons 1 and 
2. 

L.N.Lipatov [^ applied the quantum inverse scattering method to solve equations for wave 
functions of compound states of n reggeized gluons. The eigenvalue problem related to the 
operators (10) arises, because the position of the singularities in uj of the t-channel partial wave 
/jj^} is determined by their eigenvalues. 

Lipatov proposed to diagonalize the N reggeon problem by establishing the correspondence 
between the operators (p^) and the XXX Heisenberg model. He has noticed that the operator 
Hg has two equivalent representations: 

Hg = af ^ logxis^i + d:^^ logxis^s + logai^s - 2V'(1) = (11) 

Xi2 log 9i922;];2^ + 2 log Xi2 - 2^(1) 

and therefore the transposed operator can be represented in two ways as follows: 

{HgY = aiiogxi2ari + 52iogxi2a2-^ + iogaia2 - 2V'(i) = did2HG{did2)-\ (12) 

On the other hand we have 

{HgV = X12 log did2X^^ + 2 log X12 - 2V'(1) = x^^HgxI^ (13) 

Then we can deduce 

[Fg;COO] = 0, C''^ = xl^d^d2 (14) 



This equation expresses the fact of conformal invariance of Hq- Indeed, this operator coincides 
with Casimir operator of SL{2) of zero conformal weights C^^ . Then it is reasonable to denote 
Hg as H^^. 

The hamiltonian 

-f^l23 = H^2 + ^22, + -^31 = (15) 

3f ^ logX3iXi29i + d^^ loga;i2X2392 + d^^ logX23a;3i53 + 21og5i5293 - 6-0(1) = 

a;i2 log 5i52a;{2^ + X23 log 9293x^3^ + X31 log did^x^l + 2 log xi2a;23a;3i + 6-0(1) 

2 
corresponds to the exchange of 

three reggeized gluons (Fig.l). 

The transposition gives 

{Hi2'if = 5i 5253^123 (9192^3)"^ = (xi2X23a:^3i)"^-f^i23a;i2a;23a;3i, (16) 

i.e. it commutes with the operator A^: 

[Hg; A3] = 0, ^3 = xi2X23X3i9i9293 (17) 

Notice, that operator ^3 is the commutator of partial Casimir operators of the chain links: 

^3 = [C?2°, C3°?] = [C2°3°, C?2°] = [Cl'^CTs] (18) 




Therefore eq. (|17|) is a consequence of (|18|) and the Jacobi identity. It is easy to check, that the 
_/V-reggeon hamiltonian 

N 

Hn = \^Hi^ij^i, Hn,n+i = Hn,i (19) 

commutes with 

N-l N 

^^=(n^M+i)(n^j) (20) 

It can be checked also, that H^ commutes with Cat, 

Cn= E 4^*^i- (21) 

l<i<j<Af 

The eigenvalue problem: 

i^A(a:;i, 3^2, 3:3) = a^/\{xi,X2,X3), C"^ a{xi: X2, X3) = A(l - A)'I'a(xi, X2, X3), (22) 

can be considered instead of the original eigenvalue problem: 

-H'i23^A(a;i,X2,X3) = -Ei23*A(a;i,a;2,a;3), (23) 

which looks much more complicated. 

The equation (|2l|) expresses the conformal invariance of the A^ reggeon system and suggests 
the connection with the isotropic sl{2) Heisenberg model with A^ cites and cyclic boundary 
conditions, because Cn can be identified with Casimir operator of this symmetry algebra. 



Lipatov has related to each !site of the Heisenberg chain the following Lax operator 

where the SL{2) spin operators 5f 

[Sf; Sf] = ±6,,Sf, [St -,3;] = -26,,S^ (25) 

are represented as differential operators 

St = x^i^^ + 2AiXi, S^ = Xidi + Ai, S; = d^ (26) 



and A is the spectral parameter. The Lax operator (BTI) satisfies the Yang-Baxter equation. 
Therefore the trace of the monodromy matrix 



N 



T{\) = \{U{\) (27) 



is the generating function for the set of N mutually commuting differential operators Qk'- 

Jl<J2<---<Jfc 

Now it is plausible enough that all operators Qt commute with the Hamiltonian and the latter 
can be represented as a function of Q^- Explicite calculations for low N cases and additional 
mathematical arguments [^ cofirm this. Moreover, the relation between the high-energy QCD 
kernels and the XXX Heisenberg spin chains turns out to be much deeper. Namely, this 
connection can be extended to the case where fermions are incorporated. The ideas of conformal 
symmetry in Regge asymptotics |10] have been developed in application to the fermion exchange 



in |n| 



3 Review of integrable chains 

Let us review the main facts of the theory of integrable systems, concerning the XXX Heisenberg 



magnet. In our review we shall follow Sklyanin's work |12|. The phenomenon of the integrability 
of quantum systems can be understood by means of their relation to the linear ones via separation 
of variables. Namely, the quantum system is integrable, if its non-linear equations of motion can 



be represented as the zero-curvature conditions of some integrable linear system [13|. Physically 
this means, that the interaction of such systems reduces to the elastic scattering and the only 
result of it consists in the exchange of quantum numbers (momenta etc.) of the scattered 
particles. Accordingly, the 5-matrix of the theory is factorized into the product of blocks. 



corresponding to 2 ^ 2 scattering and also 1 ^ 1 in the presence of a boundary |14]. 

A set of (annihilation) operators Za{X) satisfying the Zamolodchikov algebra: has been 
proposed for an algebraic description of the factorizable scattering. 

Z,(A)Z5(m) = Sab,cd{>^ - fi)Zdifi)Z,{X) (29) 

where S is n^ x r? matrix. The consistency condition of this system, which follows from the 
associativity property of the triple product Za^{\\)Za2{X2)Za^(\'i), that is the Yang-Baxter 
equation for the 5-matrix: 

Sa.,a^ (Aj - \k)Sa.,ai (Aj - \)Sa^ai (Afc - A;) = (30) 

Sakaii^k - k)Sajaii^j - A/)S'ajafc (Aj - Afc). 



Extending this algebra by adding n conjugated (creation) operators Z|(^), one gets the Zamolod- 
chikov - Faddeev algebra: 

Za{X)Zl{fl) = 6ab6{X -fl) + Zl{fi)SacMi^ " /x)Zd(A), (31) 

or in matrix notations: 

A{X) A{fi) = Ai{\)A2{fi) = SuiX - fi)A2{fi)Ai{X), (32) 

4(A)4(m) = 4{f,)Al{X)Sl{fi - A), 
yli(A) 4(m) = /i5(A - /i) + 4{fi)Si2{fi - X)A2{X), 

where A{X) and A^{fi) are the column (Zi(A), ..., Z„(A))* and the row (Zj(;u), ..., Z^(/i) corre- 
spondingly, subscripts refer to the corresponding isotopic spaces C" (8) C" = ^1^2 and 5*21 = 
7^512^, 5*12 = VSu, 'P is the permutation operator in C" (^C". The complete scattering matrix 
'S'({Afc}) of the M particle is factorized then into the ordered product of M{M — l)/2 two-particle 
jS-matrices ([2^). For example the S'-matrix of the j'-th particle on the other M — 1 particles is 
given by t{Xj; {Am}), i-e. the particular value of the transfer matrix for A = Aj: 

t{X;{X^n^})=traT{X;{Xrn}) = tral[Sak{X-Xk), (33) 

k 

The trace in this expression is taken over the auxiliary space Va, while the transfer matrix 
acts in the quantum space CSi^^Vfc- In the framework of the quantum inverse scattering method 
(QISM) 1^, 12 1 instead of the original non- linear problem the auxiliary linear one is considered: 

— r(A, x) = L{X, x)T{X, x) (34) 

ax 

or 

T(n + l,A) = L„+i(A)r(n,A) 

in discrete case. 



This is the Lax operator of the QISM |17|. The solution of (p4): 

T{X,x) = PeMf L{X,y)dy), (35) 

T(n,A) = L„(A)L„_i(A)...Li(A) 

in discrete case, defines the monodromy matrix T{X). 

Its entries are the new variables (the quantum scattering data), which commutation relations 
are defined by 

n 
Yl ^ii2,jij2(.^ - fJ')Tj,kAX)Tj2k2{fJ') = (36) 

iij2=i 

n 
/ . ^J2J2(/^)^Jljl(A)-Rjijr-2,A:ifc2(A - /i)- 

iij2=i 

We see that integrable systems are specified by the i?-matrix, which acts on C" 0C" and satisfies 
the Yang-Baxter equation 

n 
/ . ^hi2,jlJ2\X)'njiis,kij-s[X + fJ')Rj2J3,k2k3[N = (37) 

ilj'2,J3 = l 
n 

/ J ''^i2i3,J2J3\f^)'^hJ3,jlk3\X + f^)''^jlJ2,klk2\X)- 

ilj2,J3 = l 



In general the i?-matrix depends on the spectral parameter A and other parametres. Although 
there is no complete mathematical theory of the Yang-Baxter equation, a variety of solutions 
are! known as well as different fields of their application. They are classified by the Lie algebra, 
its irreducible representations, and the spectral parameter dependence: rational, trigonometric 
and elliptic ones |^]. Given a solution -R(A) one can define the quadratic algebra Tr of n x n 
matrix elements Tij, which is generated by eq. (|36|). The associative algebra Tr realizes the 
representation space of a quantum integrable system. The commutative integrals of motion are 
t{X) = trr(A), which follows from (^), taking the trace of T1T2 = R12 T2T1R12. The algebra 
Tr posesses the comultiplication property: if Ti(A) and T2(A) are two representations of Tr in 
the quantum spaces Vi and V2, then the matrix 

Tik{\) = Ti,ij{X)T2,jk{\) (38) 

is a representation of Tr in the tensor product Vi (gi V2. This property allows to represent T(A) 
as a product of elementary representations, the so called Lax operators Li{\). It follows from 
SL{2) symmetry of the i2-matrix that an arbitrary constant d^d matrix K provides the simplest 
representation of the algebra Tr. This algebra has a central element, the quantum determinant 
ofr(A): 

A(A) = detgr(A) = D{\ + r]/2)A{\ - 77/2) - B{\ - r?/2)C(A + r//2 = (39) 

A{\ - ii/2)D{\ + 7?/2) - C{\ - 7]/2)C{X + 7?/2) = 
^(A + 7]/2)D{X - 77/2) - B{X + 77/2)C(A - r?/2) = 
D{X + r//2).4(A - r?/2) - C{X + 7]/2)B{X - r^/2) = 

which has the following remarkable properties 

detgTi(A)r2(A) = detgTi(A)detgT2(A) (40) 

and 

detgi^ = deti^. (41) 

The next representation is given by Lax operator, mentioned above, which takes the especially 
simple form for the XXX spin chain: 

L{X) = X + r^J2S^a^={^^lf ^^^„o) (42) 

where operators Sa belonging to some irreducible representaion of sl{2) have commutation 
relation (p^. Note that the i?-matrix itself can be choosen as a Lax operator, if the auxiliary 
space is two-dimensional. We have 

detgL(A) = A2-??2(c + i/4), C = {Sy -'^{S^ S' + S' S+). (43) 

Since the L(A)-operator, being the elementary representation of Tr, satisfies the Yang-Baxter 
relation and the i?-matrix depends only on the difference of the spectral parameters, the shift 
L(A) -^ L[X — uj) defines an automorphysm in Tr: 

R{X - cJi + 5f cj„) (A - u;2 + S^a/s) = (A - a;2 + S^a^)(A - wi + Sfaa)R (44) 

Separating the terms, linear in A in this equation one deduces that the i?-matrix is SL{2)- 
invariant 

[R; Sf + S^] = (45) 



and depends only on difference u;i2 = uJi —002- Tlie S'L(2)-invariance implies that the i?-niatrix 
has to have the form 

R = Y^p,{u:^2)Pj: (46) 

where Pj are the projectors corresponding to the decomposition of the tensor product of two 
initial representations into the sum of irreducible representations labelled by spin j. 

Furthermore the part of eq.(^), which contains no A gives for pj{oji2) the recurrence relation: 

Wl2 + 7?(i + l) , . , . 

P,+iK2) = ^;;^^3^^^-^/.,(^i2), (47) 

which determines i? up to a scalar factor. 

Particular solutions of the Y-B equation have properties, which are important for different 
applications, but which are not necessarily valid for a given solution: 
regularity 



P-symmetry 



T-symmetry 



unitarity 



crossing symmetry 



quasiclassical property 



i?(0) = p(0)l/2pi2 

VRi2{X)V = R2i{\) = Ri2{X) 
R%''{X) = Ru{X) 

Ri2{\) = R2i{-X) = p{X)I 

i?i2(A) = y(i)i?l|(-A-7?)F(-/ 

i?(A,r/) = / + ??r(A) + 0(r/2), 



Here the superscript t denotes matrix transposition, r(A) is the classical i?-matrix, p{X) is an 
even scalar function, r/ is the crossing parameter and V determines the crossing matrix M = 
V^V = M*. The quasiclassical property gives rise to the direct connection of the quantum model 
to the corresponding classical one. Many i?-matrices have only the combined PT-symmetry: 
-^12 (-^) ~ -R2i(A). The regularity is used to extract from t(A) the local integrals of motion. 
Thus the general solution of (E3) is given by 



T{X,u}) = KLn{X-ujn):.L2{^-uj2)Li{X-uji)=(^^^^^^^ ^J^j ) . (48) 

because any permutation of the multipliers gives the equivalent result in the algebra Tr. Notice, 
that the L-operator is noting else than a i?-matrix, acting in auxiliary and quantum spaces 
C^^Vi-. Li{X) = Rai(,X). 

The corresponding quantum determinant is 

N 

A(A) = det,r(A) = det KlliiX - coif - if{C, + 1/4)) (49) 



Now, it follows from 

i?i2(A - /i)r«(A,cJ)r(2)(^,cj) = T^^\fi,u)T^^\X,CJ)Ri2{X - /i) (50) 

that 

[t{X,u);t{X,uj)] = 0, (51) 

where t{X,Cd) = tiT{X,ij), The trace is taken over the auxiliary space. 

Among the integrals of motion ( |5l] ) we look for local ones, i.e. quantities H^^' k = 1,2, 3, .., 
which can be expressed as the sum of local operators, 

N 






^^'=^=E^J^,....-.+i (52) 



(k) 

The periodicity, A^ + 1 = 1, is supposed. The local densities Hl^_^ i-k+i should involve only 
k adjancent spins 5j, 5'j+i, ..., Si-k+i- An important case when such local integrals exist is that 
of the homogeneous spin chain, corresponding to equal spins Aj = A and zero shifts a; = 0. It 
has the important property of translational invariance. The corresponding i?-matrix is regular. 
The similarity transformation 

USfU-^ = 5f+i, USnU-^ = }CiSi}C{\ (53) 

where /C permutes the boundary matrix K and with the Lax operator Li: KLi{X) = K,^^L{X)KK,. 
This transformation generalizes the ordinary translation for the periodic chain [K = 1) to the 
twisted periodic boundary condition, specified by the matrix K and U 7^ 1 in contrast to the 
case AT = 1, when operator U takes the especially simple form: U = Vi2'P23---'Pn-in- The 
unitarity of U allows to represent it in exponential form 

U = e'^, (54) 

where operator P has the physical meaning of the total momentum of the chain. The hamiltonian 
of the model then acquires the form 

, N , 

H = ^*(A)|a=0 = E -^^i,^+lR^,i+lW\\=0■ (55) 



Faddeev and Korchemsky have shown ||9|] that the A-reggeon Hamiltonian (|15D, correspond- 
ing to the homogeneous chain can be obtained in this manner. Unfortunately, the i?-matrix, 
corresponding to the inhomogeneous chain (see below) posesses no regularity property and the 
corresponding Hamiltonian cannot be related to the derivative of the transfer-matrix in a simple 
way. 

The analysis of the integrable systems is modified upon imposing boundary conditions dif- 
ferent from the periodic ones. This case is related to the factorizable scattering of particles with 
internal degrees of freedom on a half- line [^] . The algebraic description involves additionally a 
boundary operator B into the ZF-algebra [p!g|]: 

Za{X)B = Kab{X)Zb{-X)B (56) 

Then the two particles factorizability gives rise the reflection equation (compare (p9|): 

Si2{X - fi)Ki{X)S2i{X + fi)K2{fi) = K2{fi)Si2{X + fi)Ki{X)S2i{X - ^) (57) 



in addition to the Yang-Baxter equation (pQ). The reflection matrix has the same properties as 
the it!-matrix, regularity: K{0) = /; unitarity: A'(A)Ar(— A) = /, T-symmetry: A'*(A) = Ar(A); 
the crossing symmetry is more elaborated and it involves the ^-matrix as well 111 



Then the boundary operator B can be constructed by 

B = e^p{ f ^{X)dX) (58) 

from the combination, 

<A(A) = Za{-X)Kab{X)Zk{-X), (59) 

which is a "local" field (j){X): [(l){X); (p{fj.)] = 0. Due to the 5-L(2)-symmetry of the 5-matrix, 
the corresponding i^-matrix can be transformed K ^ K' = GKG~^ with arbitrary G and the 
general solution of the reflection equation (57) for the rational case, which we are interested in, 
is 

K{X)=iI + XE, E^=I. (60) 

The reflection equation has an important covariance property: if T(X) and K(X) satisfy the 
relations (|^) and (|^) then K'{X) = T{X)K{X)T{-X)~'^ is also a solution of (|8^), provided 
the entries of K{X) and T(A) commute, [Kab{X),TcdiX)] = 0. The proof follows easily by the 
substitution of K'{X) into ( |80[ ) and by using the fundamental Y.B. relation in the different form 



r(2|(-/i)i?i2(A + A^)T(i)(A) = r(7|(A)i?i2(A + /i)r(2)(-/i). 

If the matrix T(A) is constructed as an ordered product of N independent Lax operators, 
then K'{X) can be interpreted as the monodromy matrix of N site lattice model with a boundary 
interaction described by the operator valued entries of the matrix K(X). It is called Sklyanin's 
monodromy matrix. The corresponding transfer matrix is defined as the trace 

t{X) = trK{X)T{X)K{X)T-\-X), 

where the matrix -^(A) is any solution of (p^), corresponding to the other boundary, is commu- 
tative 

[t(A),t(/i)]=0. 

In the context of the Heisenberg chain equation ( |57|) takes the form: 



i?i2(A - /.)K(-)(A)i?*^2*'(A + ^^)K^2)(^') = (61) 

A- (^)i?i2(A + /i)K(-)(A)i?if (A - f,), 
Ru{-X + /x)(i^+ )*nA)M(-Ji?*^2*'(-A -I,- 2r/)M(i)(K+))*^(/x) = (62) 

(K+))*'(/")^^(i)^i2(-A - /x - 2r/)M(-J(K+))*HA)i?if (-A + /i) 

where M is crossing matrix, defined above. In practice, if i^^(A) is a solution of (61) then 
K^{X) = {K^ {— X — ri)Y M is a solution of (62). The eq. ( p7[ ) has an important covariance prop- 
erty: if T(A,a;) and K±{X) satisfies the relations ( p7|) and (61), (62) then Sklyanin's monodromy 
matrix: 

C/(A,w) = T{X,uJ)K-{X)f{X,(:J), (63) 

where T{X,uJ) = RNa{X — CL'Ar)...i?2a(A — u;2)-Ria(A — LOi), (cr. with (^)), satisfies the relation 

RuiX- fi)U^i){X,C3)R\\'^X + fi)U^2MCo) = [/(2)(/i,t:5)i?i2(A + /x)C/(i)(A,c:;)i?*f (A - /i). (64) 
Indeed, we note that unitarity and crossing symmetry together imply the relation 

M^,)R\li-X - r,)M^^]R%{X - r,) = p{X). (65) 

Futhermore, we see that unitarity implies T(A,a;)T'(— A,u;) = np(A — ^i)- Therefore, up to a 
scalar factor, T(— A,d;) is the inverse of T(A,a;).. 



The coininutativity of the transfer matrix t(A,ti;) impUes integrabihty of the open quantum 
spin chain with the hamiltonian [^: 



whose two-site terms are given by 



Hii+i — --7r'Pii+iRii+i{X)\\=o (67) 



in the standard fashion. 



4 Closed XXX Heisenberg chain 

There are two ways of including fermions. The first corresponds to considering closed Heisenberg 
chains with different spins, i.e. contain the operators Gfg and Hp together with Hq- This 
case arises in amplitudes with the exchange of two adjoint fermions and one gluon. The same 
hamiltonian also describes the exchange of three fermions in the fundamental representation of 
5*^^(3) if not all three have the same helicity. In the first case the Regge singularity is near j = 
and in the second near j = —1/2. The seconds way of including fermions will be considered in 
the next section. 

We consider the conformally covariant operator obtained from Hp by substituting to = 0, 

H°2 = afMogxi29i + logj;i2+log 9192-2-0(1) = Xi2 log9iX];2^ + log92 + 21og j;i2 -2-0(1) (68) 
with Ai = 0, A2 = 1/2 and the conjugated operator 

H2° = 92'Mogxi292 + logj;i2+log9i 92-2-0(1) = X12 log (92X];2^ + log 9i + 2 log 2:12 -2-0(1) (69) 
with Ai = 1/2, A2 = 0. We have also Hp which should be denoted by 

H^^ = 21ogxi2 + log 9192 = xi^H^\i2 (70) 

with the weights Ai = A2 = g- We have used the identity 

(Xi292)"^ = 92'^logXi292 - logXi2 

The operator H^^ is selfconjugated. 

{H22f = H22 (71) 

and for if" 2 we have 

(H^^f = di logxi29f 1 + logxi2 + log9i92 - 20(1) = diH^^^^ = (72) 

= ^12^ log9ixi2 + log 92 + 21ogxi2 - 2V'(1) = VuxiiH^^xuVu 

Therefore, 

[H'>^;ri2Xudi]=0 (73) 

Taking into account that 

(Vuxudif = VuxudiVuxudi = X2i92xi29i = -C°^ (74) 



and comparing with the general expression 

(.A1A2 ^ (^- ^ ^- )2 ^ ^2^9^92 + 2xi2(Aia2 - Aa^i) + (Ai + A2)(l - Ai - A2) (75) 

we can conclude that indeed under confornial transformations the operator H^ transforms co- 

variantly with weights 0, g- 

Consider first the homogeneous ^ 

closed chain, which consists A 

out of three fermions. 




The corresponding hamiltonian is 



Fig.2 



111 1 11 11 



1^222 TJ2 <_ TJ2 2 _i_ IJ2 2 ClRX 

^123 --^^32+ -^21 +-^^13 - I'OJ 

2 log xi2X23a:3i + 2 log 818283. 

Conjugation arguments as described above give nothing for ^^22^ while the commutators of 
particular Casimir operators provide us with two conserved currents: 

^ii 1 = 2:12X23X3151 5253 + 1/2(X23(X31 - Xi2)5253 + (77) 

2 2 2 
X12(X23 - X3i)5i52 + X3i(xi2 - X23)535i) - l/2(xi253 + X235l + X3i52). 

is a third order differential operators and 

D^l 1 = {xl^8i82 + X^35253 + xiiOsSi - X3ia3 + X1252) + (78) 

2 2 2 

(X3i - Xi2)5i + (Xi2 - X23)52 + (X23 - X3i)53, 

is a second order one. The latter reflects the confornial symmetry of the system. The relations 

[h'4Kd['U]=0 (79) 

2 2 2 

and 

[H^iKD['L]=0 (80) 

2 2 2 

can be checked by direct calculations. These conserved currents can be obtained also in more reg- 
ular way. They appear as coefficients in front of A^ and A'^ in the monodromy matrix expansion 
i(A), where 

i(A) = tr{Ly\X)Ll^\X)Ll^\X)) = ufl , + XD^P, , + X^ + 1/4 (81) 

2 2 2 2 2 2 

and the Lax operators Lf'{X) are defined in (p^) with Aj = 1/2. 

Let us consider now 
the inhomogeneous closed chain, 
which consisits of 

one fermion and two gluons. 

The corresponding hamiltonian is 

Hi's = ^32° + h'J + 4° = (82) 

21ogXi2X23X31 + 2 log 51 52^3 + 53^^X33^ - X^^) + 82^{x^i - X^^). 




Transposition arguments do not work here. However, for closed chains with three sites there 
is always a third order differential operator, commuting with operator of total spin of chain C. 
Indeed, the result for the commutator 

j^AiA^.^AsAij ^ [xl^did2 + 2xi2(Aia2 - Aa^i); xii^sai + 2x3i(A3ai - Ai^g)] = (83) 
a;i2a;23a;3ic^i9253 + A3Xi2(x23 - x:ii)did2 + A2X3i(xi2 - X2:i)d^di + AiX23(a:3i - a:i2)9293 

-2(A2A3X239i + A3AiX3ia2 + AiA2a;i2a3) = Z^^AaAa 

(2) 

is symmetric under cyclic permutation of (123). Therefore its commutator with -D^ A A — 
0^2 ^ + ^23^ ^ + ^31^ ^ vanishes due to the Jacoby identity. However, the hamiltonian Hf2^ 

(2) ■ (3) 

commuting with D\' does not commute with D\' : 

[Hi23,D^^^]= (84) 

[(21ogXl2X23X31 + 21ogaia293 + 53^^(2:23^ - 2:31^) + d^^'^ix];^ - X^s)), 

(xi2X23a;3i3i5253 + l/2(xi2 - xl^)d2d3)] = -0^^x12X3182 - ^d2'^xl^x{2d3. 

These operators D^"^' ,D^^' also appear as coefficients in front of A^ and A*^ in the monodromy 
matrix expansion t(A), where 

t(A) = tr{L^^^{\)Ll{\)Ll{\)) = D^^'^ + AD^^^ + A=^(2A - 1). (85) 

The next inhomogeneous chain is r) 
the one with one gluon jQ 

and two fermions: 

i / 

Fig.4 
The corresponding hamiltonian commutes with the Casimir operator 

[(a:i29ixi292 + a;2353a;2392 + a;3i3i (93X31); (21ogxi2X23a;3i + 21ogai9253 + 92"^(a;f2^ - ^23^)] = 0, 

(86) 
and does not commute with the next current 

[(2;2353a;3iaixi2<92 - -^x^ndid^x^i); (21ogXi2X23a;31 + 21og9i9293 + 92"^(X];2^ - X23^)] = 

-52"^X3i(x^3^5i -Xj'2^a3)x3i. (87) 

The same happens for the longer closed chains. The hamiltonians for that chains commute 
only with the first current, the Casimir operator. This means that in order to describe the 
integrable model, the corresponding hamiltonian should be modified. Indeed, the inhomogeneous 
chains (fig. 3,4) have to be considered as the chains with impurity. Probably, the expression for 
H^2 has to be changed slightly. 

5 XXX Heisenberg chain with open boundary 

The second way of including fermions is to build an open chain with the fermions (now in the 
fundamental gauge group representation) at the ends. This open chain corresponds to the Regge 
exchange with meson quantum numbers in the i-channel. 




Let us consider the hamiltonian ^ „ o corre- 

• — o — o 

Fig.5 
spending to the following chain 



'l23 ~ -"32 "t" -"21 



2 



log a;i2 + ^2 ^ log X32a;2i92 + log 818283 + ^3 ^ log X3293 = 

log 81 + X32 log 8283X^2 + ^21 log 52X^/ + 2 log X21X32, 

The transposed Hamiltonian has form 

(-f^l23 )^ = 8283111238283 = 7'i23(a;322;2l)"^-H'l232;322;21^123- 

Then H commutes with 

-00 
^123 = ^123a:^322;21<92(93^123- (89) 

The permutation operator P123 maps the sites 1,2,3 of the chain into the sites 3,2,1 correspond- 
ingly. 

In a similar way the conserved operator of the highest order in the derivatives can be obtained 
for the open chain with N — 2 gluonic operators and fermions of opposite spin at the ends: 

Af-l 



Hlt.'.N = 92^ log ^21^2 + log X21 + log 92^1 + J] H^^_^,, (90) 



tO...O 

= ^2 " log 2:21 0*2 + iog a;2i + log ^201 + 

4 = 2 

commutes with the charge operator 

A12...N =Vi2...NXNN-1--X2i8n8n~1--82. (91) 

For simplicity we put the reflection matrices K equal to unity. From eq. (^6|) one can see 
that this implies the transfer matrix to be an even function of the spectral parameter. In our 
case we obtain up to insignificant c-number terms 

t{X) = (4A2 - 1)(I)W - .2i^Po - (A^ - o^Di^Po + (^' - -^2)')' (92) 



where D\l^ and -DW are the same as in (|38D for the closed chain and 

-^iOD ^ (^122^23'9l92 + 2;i2X23(xi2 - X23)5i + Xi2xl.^82 + X23(a;i2 " 2:23)9253. 

This coincides with the square of the operator ([75|), obtained by transposition arguments. 

We put here wf = wl = i^i + 3/4 in order to achieve that the second order operator D^"^' 
coincides with C*^^". Now the operator D^"^' commutes with hamiltonian. Thus, for the open 
chain (Fig. 5), described by the hamiltonian H123 the eigenvalue problem like eq. (|6^) can 
be replaced by the one with the same Casimir operator and with the fourth order differential 
operator D^'^' instead of D^^' = A as for the closed chain. 

So we see that for open chains we are able to find conserved currents and its number is 
just enough to solve the! eigenvalue problem for QCD reggeon interactions. Moreover, for 



the open chain Fig. 5 with one ferniion inside again everything is all right. The haniiltonian 
o -, r, and the transfer matrix are given by 

o — • — o 

Fig. 6 

hS? = hI^ + H^^ = 2 log ^12X31 + 2 log dld2d^ + d^' log Xi2 - Sg-l log X31 (93) 

t(A) = 4A(A + 1)X31 91X2291X319293 + 4A2(A + 1)'(X239293 + X3l9iZ3l93 + Xi29iXi292 (94) 

The case of the open chain with two fermions and one gluon can be considered as well. 
For the chain „ -. „ 

• — o — • 



Fig. 7 
we have two conserved currents: 



[xi2X239l 9293X12X2392; (21ogXi2X23 + log9i9|93 + 82 ^(Xi2^ - Xgg^)] = (95) 

[(xi29ixi292 + X2393X2392 + X3i9i 93x31); (21ogxi2X23 + log9i9|93 + d2^{x{i - X2i))] = (^96) 



Another chain with two fermions and one gluon is: „ 



o 



Fig.8 

We have following relations for this chain 

[xi292X2393X2392Xi29i; (2 log X12X23 + log9i9|93 - d{'^x{2)] = (97) 

[(xi292Xi29i + X239293X23 + X3i93X3i9i); (21ogxi2X23 + log9i9|93 - d{^{x{2)] = (98) 

For completeness 
let us consider also 

the open chain with three fermions o i o 



Fig.9 



we can deduce following relations: 



[(9lXl292X2393X2392Xl2 + Xl292X2393X2392Xl29l); (21ogXl2X23 + log 9l 9^93)] = (99) 

[(xi29i 92X12 + X239392X23 + a;3i9i 93x31); (21ogxi2X23 + log9i9|93)] = (100) 

Thus, all possible open chains with three sites corresponding to gluonic or fermionic reggeons 
are integrable. The conserved currents can be derived from the transfer matrix with boundaries 
represented by unity reflection matrices. The hamiltonians are constructed from the QCD 
kernels given above and commute with these two conserved currents given by the second and 
fourth order differential operators respectively. We believe that integrability holds for the longer 
open chains too. 



6 Conclusions 

The examples considered above, except for the difficulties with the inhomogeneous closed chains, 
show the validity of Lipatov's conjecture about the deep connection between the kernels of high- 
energy QCD scattering amplitudes and the exactly solvable two-dimensional models in the cases 
with fermions as well. In the pure gluonic this connction context has already been used for 
the solution of the Odderon problem [^^ . The case of open chain differs from the closed chain 
case by the higher order of non-trivial conserved charge; the fourth order appears instead of the 
third. 

Recently it has been shown that the considered integrable structures appear also in hard 
(exclusive and deep-inelastic) scattering ||2^, corresponding to a different limiting region of high 
energy scattering amplitudes in QCD. 
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